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The source content of low galactic latitude XMM-Newton surveys 
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Abstract. We present results from a project conducted by the Survey Science Center of the XMM-Newton satellite and 
aiming at the identification and characterisation of serendipitous EPIC sources at low galactic latitudes. Deep multi-colour 
optical imaging and spectroscopic observations have been obtained in the framework of several observing campaigns carried 
out at ING, CFHT and ESO. These observations have lead to a number of optical identifications, mostly with active stars. 
We describe the identified source content at low galactic latitudes and compare stellar populations properties at low and high 
galactic latitudes with those expected from stellar X-ray count models. 
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1. Introduction 

Thanks to their excellent sensitivity, good image quality over 
a wide field of view (30') and wide energy range (0.2-12 
keV) the EPIC cameras on board XMM-Newton allow X- 
ray surveys with an unprecedented combination of large area 
and depth. This potential was recognised by ESA in setting 
up a dedicated XMM-Newton Survey Science Center (SSC) 
to facilitate the exploitation of the XMM-Newton Serendipi- 
tous Sky Survey by providing a public archive of data prod- 
ucts and carrying out a carefully coordinated follow-up pro- 
gramme to characterise the overall X-ray source population 
(Watson et al. 2001). 

In this paper we present results from an ongoing SSC op- 
tical campaign whose aim is to provide a reference frame 
for the characterisation of all low galactic latitude EPIC 
serendipitous detections. Beyond proper identification and 
classification, our objective is to study source properties, 
quantify their relative populations and contributions to the 
overall soft and hard X-ray galactic emission. XMM-Newton 
measurements are not strongly biased by interstellar extinc- 
tion, as were earlier studies based on ROSAT observations 
restricted to the softer 0.5-2 keV band. Although similar 
in X-ray flux to the deepest ROSAT pointings, the galactic 
landscape unveiled by XMM-Newton differs from that seen 
by ROSAT. Opening new parameter spaces the present sur- 
vey offers the possibility to reveal completely new classes of 
sources. 

Many individual source types are known to contribute 
to the galactic X-ray source demography in the medium to 
low X-ray luminosity regime (10 28-36 erg s _1 ). These in- 
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elude, late and early-type stars, cataclysmic variables (CVs), 
RS CVns and various species of X-ray binaries. Other 
kinds of rare emitters such as isolated neutron stars or low- 
luminosity but long-lived evolutionary states of classical 
XRBs may be found. 

2. Observational material 

Target fields (see Table Q were selected for their good X- 
ray quality (depth and low background) and chosen in di- 
rections void of extended diffuse emission or bright sources. 
We also avoided atypical regions such as known star forming 
complexes. Mainly based on EPIC pn data, source detections 
were visually screened and positions were in most cases cor- 
rected for residual attitude errors by cross-correlating them 
with large optical catalogues. In order not to exclude extreme 
sources we retained detections in all energy bands. 

Optical data were collected in the framework of the AXIS 
project (Barcons et al. 2002) performed at the Observatorio 
del Roque de los Muchachos in the Canary Islands and from 
programmes conducted at CFHT and at ESO. The wide field 
camera at INT, the CFH12K at CFHT and the WFI at ESO- 
MPG 2.2m provided wide-field imaging in the g', r', i', z, 
R, I and Ha / Ha continuum from which high spatial res- 
olution mapping of the error circles and candidate priorities 
for spectroscopic observations were derived. Medium to low 
resolution spectroscopy was obtained on 4m-class telescopes 
using multi-fiber instruments (WHT/WYFFOS) or long slit 
spectrographs (WHT/ISIS, ESO 3.6m/EFOSC2). 

The total area surveyed so far is ~ 0.8 deg 2 . With expo- 
sure times in the range of 10 to 30 ksec we reach flux lim- 
its of - 4 10 ~ 15 erg cm" 2 s _1 in the 0.5 to 2.0 keV soft 
band and of 



3 10 ~ 14 erg cm" 2 s" 1 in the 2 to 10 keV 
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Table 1. Galactic target fields 



Field 


RA 


Dec 


1 


b 


Gal Nh 


WR 110 


18h07 


-19°23' 


10.8° 


+0.4° 


5.5 IQr 1 ' 2 


Ridge3 


18h27 


-11°29' 


20.3° 


+0.0° 


1.1 10 23 


G21.5-09 


18h33 


-10° 34' 


21.5° 


-0.9° 


9.0 10 22 


Z And 


23h33 


+48° 49' 


109° 


-12.1° 


1.4 10 21 



hard energy range. The present survey can thus be qualified 
as medium sensitivity in being about 10 times shallower than 
the Chandra deep galactic survey of Ebisawa et al. (2002) but 
10 times deeper than the medium ROSAT survey of Morley 
et al. (2001) or than the faint ASCA galactic plane survey of 
Sugizakiet al. (2001). 

In order to compare stellar population properties in the 
galactic plane with those observed at high galactic latitudes 
we also use stellar identifications obtained in the course of the 
high \b\ AXIS programme (Barcons et al. 2002). The sample 
consists of 34 stars extracted from 18 XMM-Newton fields 
(|6| = 23°— 80°, average value 50°) in a total survey area of 
2.3 deg 2 . The flux limit is 2 10 " 14 erg cm" 2 s" 1 in the 0.5 
to 4.5 keV band slightly above that used in the galactic plane. 



3. Identification Strategy 

Identification of XMM-Newton sources in the Milky Way can 
be made difficult by optical crowding. Matching the pattern 
of EPIC sources with the USNO A-2 catalogue and then cor- 
recting for residual attitude errors can shrink the 90% con- 
fidence level error radius down to ~ 2.5". However if field 
stellar density is too high, no such adjustment can be made 
and the 90% confidence level error radius remains of the or- 
der of 4.5". Based on our wide field imaging we estimate that 
at b ~ 0° the probability to find at random a star brighter than 
R~ 16 - 17 in a typical EPIC error circle is low enough (5%) 
that we can accept the identification on the basis of positional 
coincidence. However, in many cases this argument cannot 
be used and we have to collect optical spectra of typically 3- 
4 candidates brighter than R ~ 22 before possibly finding a 
reliable identification. At the sensitivity of our survey, stellar 
counterparts fainter than R ^ 16 - 17 are mostly late K-M 
stars. Balmer emission is found in all active M stars and in 
a fraction of the most active K stars. Cataclysmic variables 
and Be/X-ray binaries are also Ha emitters. Our narrow band 
Ha and Ha continuum imaging together with broad band 
colour information can thus be used to efficiently prioritize 
optical candidates for spectroscopy (Herent et al. 2002). We 
also discarded sources with hardness ratios consistent with 
those expected from an absorbed background AGN and with 
only faint optical candidates in the error circle. This selec- 
tion is efficiently done in directions of high absorption and 
still leaves a large range of possible Nh (or intrinsic spectral 
hardness) for galactic objects. 
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Fig. 1. Optical spectrum of a faint accreting candidate in the 
field of WR 110 

4. Optical identifications 

We list in Table |2] the statistics of optical identifications in 
each field and for the entire galactic sample. The limiting 
magnitude defined here as the faintest magnitude at which 
we can detect emission lines typical for a CV or a Me star 
depends on the instrumental setting and is R ~ 19 for the 
G2 1.5-09 and Z And fields and R - 21 for the Ridge 3 and 
WR 110 fields. The vast majority of optical identifications 
are with active coronae. Because of the relatively high galac- 
tic latitude and correspondingly low Nh it was still possible 
to spectroscopically identify two AGN in the Z And field. At 
lower latitudes, galactic absorption becomes too large for our 
instrumentation. The only low b extragalactic source is in the 
field of G21 .5-09 (XMMU J183225 .4-103645, Nevalainen et 
al. 2001) and was identified as an X-ray bright cluster of 
galaxies seen through over 50 magnitudes of absorption in 
the visual on the basis of X-ray source extent and energy dis- 
tribution. 

4.1. Accreting candidates 

In three cases, we detect X-ray emission from a Be star at a 
level which is slightly above that expected for 'normal' early 
type stars (Cassinelli et al. 1994). In addition to the case of 
SS 397 already mentioned in Motch (2000), we find excess 
X-ray emission from USNO 0750-13549725 which is the 
brightest star in the open cluster NGC 6649. At a distance 
of 1.6 kpc and with E(B-V) = 1.38 (Turner 1981) the X-ray 
luminosity is ~ 4 10 32 erg s _1 (0.5 - 7.5 keV). Finally, we 
detect Ha emission from a faint red star (R = 22.0, R-I=2.6) 
in the field of WR 110. The absence of TiO features and the 
colour index indicate high interstellar absorption (see Fig.[0. 
If the counterpart is an intrinsically blue object then Av is of 
the order of 12.5 orN H = 2.2 1 22 cm" 2 , i.e. about half of the 
total expected galactic absorption in that direction. The high 
interstellar reddening and accordingly large distance rules out 
a CV nature. A Be star identification seems likely and the star 
would then be at ~ 12 kpc implying an X-ray luminosity of 
1.3 1 33 erg s" 1 (0.5 - 7.5 keV). Such modest X-ray lumi- 
nosity excesses may be in part due to errors on the reddening. 
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Table 2. Statistics of optical identifications at low galactic latitudes 



Field 


G21.5-09 


Ridge3 


WR 110 


Z And 


Total 


b 


-0.9° 


0.0° 


+0.4° 


-12.1° 




Area 


0.27 


0.18 


0.18 


0.18 


0.81 


Number of sources 


77 


30 


58 


38 


203 


Stellar coronae 


15 (19%) 


14 (47%) 


29 (50%) 


9 (24%) 


67 


Accreting candidates 


2 





1 





3 


Extragalactic sources 


1 








2 


3 


Unidentified 


59 (77%) 


16 (53%) 


28 (48%) 


27 (71%) 


130 



We note however, that the numerous Be + WD systems pre- 
dicted by binary evolution theory could well be emitting in 
this range of luminosities (see Raguzova2001 and references 
therein). 



could be due to close binaries such as RS CVns for which 
synchronized rotation maintains a strong X-ray luminosity 
over a long time. Alternatively, a fraction of these high Lx 
stars could be misidenfifications. 



4.2. The stellar population 

At very low galactic latitudes, the fraction of sources identi- 
fied with active coronae varies very significantly with point- 
ing direction. Whereas in the WR 110 and Ridge 3 fields stel- 
lar coronae account for about half of the X-ray sources, the 
fraction is only 19% in the field of G21.5-09 in spite of the 
fact that source density is similar to that of WR 110. Such a 
variance is not unexpected and could reflect the presence of a 
higher 'local' absorption in the field of G21.5-09 preventing 
spectroscopic identification of the globally optically dimmer 
stellar population. 

Spectral types, distances and luminosities were derived 
for most stellar identifications in the WR 110 and Ridge 3 
fields as well as for the high galactic latitude sample (see Ta- 
ble|3}. 

In both 6^0° fields, the mean B- V excess of X-ray active 
stars is about 0.6 (or Nh ~ 3.5 10 21 cm -2 ). We detect active 
coronae up to 500 pc in Ridge 3 and 1 kpc in WR 110. 

At the flux level considered here (~ 4 10 ~ 15 erg cm -2 
s _1 in the 0.5 to 2.0 keV band), the low \b\ active star pop- 
ulation displays a distribution in spectral types which is not 
statistically different from that seen at high The low lati- 
tude ROSAT survey (Motch et al. 1997) also exhibits a sim- 
ilar spectral type distribution in spite of a factor 25 lower 
sensitivity (see Table [3}- At low and high latitudes, the ob- 
served distributions match well those expected from X-ray 
count models (Guillout et al. 1996, Guillout & Motch 2002). 
In particular, we are not yet in a flux regime faint enough for 
M stars to dominate at high latitudes. Deeper surveys could 
however reveal this effect. 

The distribution in X-ray luminosity of the two b ~ 0° 
stellar samples and of the |6| ~ 50° sample is shown in Fig. 
12 In the galactic plane, Lx peaks between 10 29 and 10 30 erg 
s _1 , typical for stars with ages of the order of that of the 
Pleiades (< 10 8 yr, see e.g. Micela et al. 1990). In contrast, 
the Lx distribution at high \b\ shows two peaks, one between 
10 28 and 10 29 erg s _1 and one at very high Lx between 10 30 
and 10 31 erg s _1 . The low luminosity component has a so- 
lar like X-ray luminosity typical for the relatively old stel- 
lar population preferentially detected at high latitudes. The 
origin of the high luminosity component is not yet clear. It 



Table 3. Distribution in spectral types of active coronae (per- 
cents) 

XMM-Newton : Observed and Model prediction 







M 


K 


G 


F+A 


\b\' 


-0° 


28.6 ± 6.4 


32.7 ± 6.7 


22.4 ± 6.0 


16.3 ± 5.3 






24 


25 


31 


20 


\b\ ' 


-50° 


29.4 ± 7.8 


17.6 ±6.5 


23.5 ± 7.3 


29.4 ± 7.8 






27 


19 


29 


25 


ROSAT all-sky survey 






M 


K 


G 


F+A 


\b\. 


-0° 


19.0 ± 6.0 


23.8 ±6.5 


26.2 ± 6.8 


30.5 ± 7.1 



We also show in Fig.[3]the log N(>S)-log S curve for stars 
at low and high galactic latitudes. At EPIC pn count rates of 
a few 10~ 3 cts/s (0.5 - 2.0 keV), the density of active stars is 
about 10 times larger in the galactic plane (here only Ridge 
3 and WR 110 fields) than at high b. Predictions of the stel- 
lar population models for b = 0° and 50° fit rather well the 
observed curves. One should however keep in mind the large 
density variance observed at low latitudes. As argued above, 
part of this scatter could be due to complex absorption struc- 
tures which are not yet taken into account in the X-ray count 
model. X-ray luminous old binaries may also account for a 
sizeable fraction of the population. Clearly more stellar iden- 
tifications are needed in order to build really representative 
samples. 

5. Conclusions 

Not unexpectedly, our observations confirm the well known 
concentration of young active stars in the galactic plane. Be- 
cause of the steep decline of coronal activity with age, X- 
rays highlight the young stellar population. Whereas at op- 
tical wavelength stars younger than 1 Gyr make up only ~ 
10% of the total stellar population and do not exhibit marked 
signatures, they account for about half of the X-ray detec- 
tions, the exact fraction depending on galactic latitude and 
limiting flux. The large sensitivity of XMM-Newton allows 
to sample young stellar populations up to 1 kpc or more as 
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Fig. 2. Distribution in X-ray luminosity of active coronae. Left: |6| ~ 0°. Right: |6| ~ 50° 



demonstrated in the WR1 10 field. At these distances, a num- 
ber of effects expected from X-ray stellar populations models 
become readily visible and can be used to put very significant 
constraints on the dynamical properties of young stars, espe- 
cially the evolution of scale height with age, and on the SFR 
during the last ~ 2 Gyr. Among these effects are the sharp 
decrease of number count with increasing galactic latitudes, 
the rising fraction of old stars at high b and the increasing 
fraction of M stars at faint fluxes. 
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Fig. 3. X-ray stellar logN(>S)-logS curves at low and high 
galactic latitudes. At \b\ = 0° the green (or light) plot repre- 
sents the log N(>S)-log S curve in Ridge 3 and the black one 
that in WR 110. Continuous (red) lines show model predic- 
tions for \b\ =0°and \b\ =50° 
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